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Abstract. Diagnostic length differences in a PCR am¬ 
plified fragment of the gene for byssal adhesive protein 
were used to study the zoogeographic distribution of Myti¬ 
lus galloprovincialis and M. trossulus along the west 
coast of North America and in Japan. The distributions 
of M. galloprovincialis and A7. trossulus are patchy, al¬ 
though an overall geographic pattern emerges. Al. gallo¬ 
provincialis was the only species found on either Kyushu 
or Honshu, and it was the most abundant mussel from 
Tomales Bay to San Diego, California. M. trossulus was 
the only bay mussel found on Hokkaido and in Alaska, 
and it was by far the most abundant mussel along the 
coasts of Washington and Oregon. 

Mytilus galloprovincialis and M. trossulus are sympat- 
ric and hybridize near Whidbey Island. Washington, in 
San Francisco Bay, and in San Diego Bay. A second 
diagnostic anonymous nuclear PCR marker was used to 
examine the extent of hybridization at Palo Alto. Califor¬ 
nia. At this site, genotypes appeared to be a mixture of 
M. galloprovincialis, F] hybrids between M. galloprovin¬ 
cialis and Ah trossulus, and backcrosses between the F/s 
and M. galloprovincialis. 

The discontinuity between the zoogeographic distribu¬ 
tions of these two species at about 40°-4PN latitude in 
both the eastern and western Pacific suggests that temper¬ 
ature is a factor in determining their present distribution 
and limiting their dispersal to other regions. 

Introduction 

Biochemical and molecular systematics have dramati¬ 
cally altered our understanding of the ecology, evolution. 
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and biogeography of bay or blue mussels of the genus 
Mytilus (McDonald and Koehn, 1988; Koehn, 1991; Seed, 
1992; Rawson and Hilbish, 1995). The taxonomy of the 
Mytilus edulis complex has been confused (see Seed, 
1992, for a full discussion of species definitions in this 
genus) as a result of considerable morphological similar¬ 
ity among species. For example, each species varies con¬ 
siderably in size and shape both within and among popu¬ 
lations (Seed, 1968). Today, multivariate analyses of mor¬ 
phological shape (McDonald et ah, 1991) and analyses 
of molecular genetic variation (see below) have led to the 
recognition of three sibling species: M. edulis Linnaeus 
1758, M. galloprovincialis Lamarck 1819, and M. trossu¬ 
lus Gould 1850. Early work suggested a simple geo¬ 
graphic distribution in the northern hemisphere, with M. 
edulis being found in the temperate regions of the Atlantic 
and Pacific Oceans, and M. galloprovincialis restricted to 
the Mediterranean Sea. General recognition of M. tf ossu- 
lus resulted from surveys of allozyme variation in popula¬ 
tions at higher latitudes in both the Atlantic and Pacific 
Oceans (Koehn et ah, 1984; McDonald and Koehn, 1988; 
Varvio et ah, 1988). 

The discovery of hybridization among members of the 
Mytilus edulis complex blurred the tidy distribution pat¬ 
terns mentioned above. The first indication that the distri¬ 
butions were more complex than initially proposed was 
provided by electrophoretic surveys of genetic variation 
of mussels in the British Isles (Skibinski and Beardmore, 
1979; Skibinski, 1983; Koehn, 1991, p. 134), which re¬ 
vealed the presence of M. galloprovincialis, and hybrid¬ 
ization between M. edulis and M. galloprovincialis (Ski¬ 
binski and Beardmore, 1979; Gosling and Wilkins, 1981; 
Skibinski, 1983; Skibinski etah, 1978; 1980; 1983). Stud- 
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ies also revealed that the bay mussels on the Pacific coast 
of North America were not M. edulis, nor were they a 
single species. McDonald and Koehn (1988), using allo- 
zyme markers, identified Al. galloprovincialis in southern 
California, Al. trossulus in northern California and further 
north, with a hybrid zone in and around San Francisco 
Bay. This pattern has been verified with studies using 
markers from mtDNA (Geller, 1994; Geller et al., 1994; 
Rawson and Hilbish, 1995). 

Complexity and patchiness in the distribution patterns 
of Mytilus species have been found not only on a broad 
geographic scale but also on a microgeographic scale. 
For example, ecological genetic studies have revealed 
variation associated with environmental conditions. In the 
British Isles, the abundances of A/. edulis and Al. gallo¬ 
provincialis vary with both salinity and wave exposure: 
Al. galloprovincialis is relatively more common in areas 
with high salinity and high wave exposure (Skibinski et 
al., 1983). In a hybrid zone at the salinity gradient be¬ 
tween the North Sea and the Baltic Sea, M. edulis alleles 
give way to A/, trossulus alleles as salinity decreases (Vai- 
nola and Hvilsom, 1991). 

In the past, taxonomic work on Mytilus was made dif¬ 
ficult by the lack of an easy and complete means of identi¬ 
fying species. Although identification of A/, edulis, M. 
galloprovincialis, and A/, trossulus is possible with elec¬ 
trophoretic surveys of proteins (Koehn et al., 1984; Grant 
and Cherry, 1985; McDonald and Koehn, 1988; Vainola 
and Hvilsom, 1991), it is now easier to identify species 
and their hybrids with DNA markers (Geller and Powers, 
1994; Geller et al., 1994; lnoue et al., 1995; Rawson and 
Hilbish, 1995). For example, diagnostic mtDNA markers 
were used by Rawson and Hilbish (1995) to describe the 
distributions of A/, galloprovincialis and M. trossulus 
from San Diego to Seattle. Identification of species with 
mtDNA is, however, more difficult in mussels than in 
other groups of species, for salt (Mytilid) and freshwater 
(Unionid) species have both maternally and paternally 
inherited mtDNA systems (Skibinski et al., 1994; Zouros 
et al., 1992; 1994; Geller. 1996; Hoeh et al., 1996; Liu 
et al., 1996). The difficulties of using mtDNA as a diag¬ 
nostic marker have been overcome by the use of a nuclear, 
codominant diagnostic marker, lnoue et al. (1995) de¬ 
scribed a diagnostic length difference in a gene that codes 
an adhesive protein found in byssal threads. The fragment 
may be amplified readily from genomic DNA. and the 
three species have distinctly different fragment sizes, 
making it simple to identify all three species and their 
hybrids, in which two fragments are coamplified. 

The primary objectives of this study are to describe the 
macrogeographic and microhabitat distribution of three 
species (A/, edulis, M. trossulus, and M. galloprovincialis) 
and the occurrence of their hybrids. Our sampling is more 
extensive than in previous studies, with multiple sites at 
some localities and a greater geographic extent of sample 


sites. Finally, we have used a codominant, nuclear marker 
to identify the three species and their hybrids, and em¬ 
ployed a second diagnostic marker to test for the presence 
of backcross genotypes within a hybrid zone. 

Materials and Methods 

Sample sites 

Collection sites are shown in Figure 1 and described 
in Table 1. Mussels from Tillamook Bay, Coos Bay 
(1991), San Francisco Bay, Elkhorn Slough, Morro Bay, 
and San Diego Bay were the same as those reported upon 
in Geller et al. (1994). Additional mussels were obtained 
from Whidbey Island in the Puget Sound and shipped live 
to Wilmington, North Carolina, whereas mussels from 
Yokohama, Japan, were preserved in ethanol. DNA ex¬ 
tractions and amplifications of these samples were con¬ 
ducted by JBG. The remaining samples were handled by 
BRK and JBM. Samples were shipped, either alive or 
frozen, to Boulder, Colorado, and held in an ultracold 
freezer (-S0°C) until DNA was extracted. 

DNA markers 

DNA was extracted from the mantle tissue of each 
mussel by following the protocol of Geller et al. (1994). 

In the study by lnoue et al. (1995) Me-15 and Me-16 
PCR primers for a portion of an adhesive protein gene 
produced species-specific length variants of 180, 168, and 
126 bp for Mytilus edulis, Al. trossulus, and Al. galloprov¬ 
incialis respectively. Amplifications by BRK and JBM 
were conducted in a total volume of 25 p\ using 10 mA/ 
Tris (pH = 8.3), 50 mA/ KC1, 0.01% gelatin, 200 pAl 
dNTPs, 4 mA/ MgCL, 1.5 units Taq DNA polymerase, 
2 pmol of each primer. 0.75 p\ DNA template, and auto¬ 
claved ultrapure water. Cycling conditions consisted of 
an initial denaturing step of 94°C for 45 s. This was fol¬ 
lowed by 35 cycles of 94°C for 45 s, 56°C for 30 s, 70°C 
for 90 s and concluded with a final extension step of 70°C 
for 3 min. Amplifications conducted by JBG used 10 mA/ 
Tris (pH = 8.3), 1.5 mA/ MgCL, 50 mA/ KC1, 0.01% 
Triton-X 100. 0.01% gelatin, 0.01% NP-40, 200 pM 
dNTP s, 1.5 mA/ MgCL, 0.3 units Taq DNA polymerase, 
25 pmol of each primer, 0.5 p\ DNA template, and auto¬ 
claved ultrapure water in a volume of 25 p\. Cycling pa¬ 
rameters were 30 s at 98°C followed by 30 cycles of 
10 s at 95°C, 30 s at 54°C, and 30 s at 72°C. For all 
amplifications, 10 p\ of the PCR product were run on 3% 
agarose gels, stained with ethidium bromide, and visual¬ 
ized under UV light. 

A pair of primers for a single copy anonymous nuclear 
DNA locus (Cmg-93) in a deep-sea hydrothermal vent 
vesicomyid clam, Calyptogenea magnifica (Karl et al., 

1996), was found to amplify a single band of about 750 bp 
in Al. trossulus but produced no product in either Al. 
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Figure t. Mytilus collection sites i 

edulis or M. galloprovincialis. This marker, in conjunc¬ 
tion with the diagnostic marker of Inoue et al. (1995), 
permitted a more thorough analysis of genotypes in a 
zone of hybridization between Al. trossulus and Al. gal¬ 
loprovincialis at Palo Alto, California. F, hybrids between 
Al. galloprovincialis and Al. trossulus would have two 
bands at the Inoue marker, and a product for the Cmg- 
93 locus. An individual with two bands at the Inoue 
marker and no Cmg-93 product is one of several geno¬ 
types produced by backcrossing an F, hybrid with a speci¬ 
men of A/, galloprovincialis. The presence of this geno¬ 
type clearly indicates that hybridization has extended be¬ 
yond the F, level. 

Amplifications of Cmg-93 were conducted by BRK 
and JBM in a total volume of 25 p\ using 10 mA/ Tris 
(pH = 8.3), 0.01% gelatin, 50 mA/ KC1, 300 pAl dNTPs, 
2.5 mA/ MgCF, 1.5 units Taq polymerase, 6.25 pmol each 
of primer, bovine serum albumin (0.1 pg/p 1), 0.75 p\ 
DNA template, and autoclaved ultrapure water. Cycling 
conditions consisted of an initial denaturing step of 94°C 
for 60 s. This was followed by 35 cycles of 94°C for 60 s, 
50°C for 60 s, 70°C for 60 s and concluded with a 70°C 
2-rnin extension step. Ten microliters of the PCR product 
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Japan (A) and North America (B). 


were run on 1% agarose gels, stained with ethidium bro¬ 
mide, and visualized under UV light. 

Annual and monthly (January and July) means for tem¬ 
perature and salinity for each site were taken from the 
NOAA World Ocean Atlas (1994). Data were obtained 
from the nearest point for each site from a 1° X 1° grid 
of the North Pacific Ocean. In addition, some indication 
of the degree of wave exposure was recorded for most 
of the collection sites. 

Results 

On both sides of the North Pacific Ocean. M. gallopro¬ 
vincialis predominates in the lower latitudes, and M. tros¬ 
sulus predominates in the higher latitudes (Table 1, Fig. 
1). In Japan, we found Al. galloprovincialis alone at all 
sites on the southernmost island. Kyushu, and on Honshu, 
whereas only A/, trossulus was found on the northern 
island of Hokkaido. M. galloprovincialis is most abundant 
in southern California, but only A/, trossulus is found in 
Alaska. We did not detect M. edulis in any of our popula¬ 
tion samples. 

Hybridization was not detected in the samples from 
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Japan, but hybrids were common in San Diego Bay and 
San Francisco Bay. Both M. galloprovincialis and A/. 
trossulus and their hybrids were detected at Whidbey Is¬ 
land, near a mussel farm that raises M. galloprovincialis. 
A single hybrid was found in a sample from Tatoosh 
Island, Washington. 

The patchy distribution of species is apparent in and 
around the bays of California. For instance, only M. gal¬ 
loprovincialis was present on the outer coast at Scripps 
Pier in San Diego, but just a few kilometers away, A/. 
galloprovincialis, M. trossulus and their hybrids were 
found in San Diego Bay. In San Francisco Bay, near the 
Golden Gate Bridge at North Beach, both species and 
hybrids were common (34% M. galloprovincialis, 11 % 
M. trossulus, 55% hybrids). However, deeper in the bay, 
at Palo Alto, M. trossulus was apparently absent, but A7. 
galloprovincialis and hybrids were common. The individ¬ 
uals collected from Palo Alto were analyzed with the two 
diagnostic markers to examine the extent of hybridization 
(Table 1). In both 1990 and 1996, we identified genotypes 
consistent with F, hybrids (14% and 20%. respectively) 
and backcrosses (14% and 30%, respectively) of F, hy¬ 
brids to M. galloprovincialis. 

In Isthmus Slough and South Slough within Coos Bay, 
only M. trossulus was found in 1991. However, a sample 
taken from Isthmus Slough in 1996 contained both M. 
galloprovincialis and M. trossulus. It is not known 
whether the difference between the two samples taken 
from Isthmus Slough reveals a change over time or pre¬ 
existing microgeographic variation. 

Although our survey identified microgeographic varia¬ 
tion in species distributions, an environmental cause was 
not apparent. For our collection sites, neither wave expo¬ 
sure (Suchanek, 1978) nor salinity supplied any predictive 
power to explain the distribution of species. 

Despite the lack of any explanative value in environ¬ 
mental conditions at a microgeographic scale, on a 
broader geographic scale the transitions between M. gal¬ 
loprovincialis and M. trossulus occur between 40° and 
41° north latitude on both sides of the Pacific. Data in 
Table 1 allow a comparison of temperature regimes across 
the transition zones. At the break on the eastern side 
of the Pacific, there are only slight differences in mean 
temperatures (compare sites 22 and 23). The only appre¬ 
ciable difference is a drop of about 2°C in the mean 
January temperature at the northern site (site 22). How¬ 
ever, a comparison of sites 8 and 9 in Japan reveals sub¬ 
stantial differences in temperature across the transition 
zone. The northern site not only has a lower mean January 
temperature (by about 5.5°C), but also lower mean annual 
and mean July temperatures (by 3°C and 7°C, respec¬ 
tively). A comparison of the eastern and western sites at 
the southern boundary of pure M. trossulus (sites 9 and 
22) reveals that mean January temperatures are quite dif¬ 
ferent (2.21° and 10.5°C, respectively), but that annual 


mean temperatures (11.72° and 13.74°C, respectively) and 
mean July temperatures (13.39° and I3.74°C, respec¬ 
tively) are similar. 

Discussion 

In comparison to previous studies, our finer scale geo¬ 
graphic sampling of Mytilus in California revealed a more 
patchy distribution of pure species and a greater geo¬ 
graphical extent of hybrids. Whereas previous studies had 
suggested that pure M. galloprovincialis extended as far 
north as Point Conception, we found a population of pure 
M. galloprovincialis in Tomales Bay, 700 km further 
north. Furthermore, our studies of populations on both 
sides of the Pacific revealed that the transition zones are 
at similar latitudes, and at similar temperatures, sug¬ 
gesting that temperature influences the distributions of 
these species. 

Hybridization and patchiness 

The authors of a previous study (Rawson and Hilbish, 
1995) concluded that hybridization of A/, galloprovin¬ 
cialis and A/, trossulus was rare in southern California. 
They sampled from the pier at the Scripps Institution of 
Oceanography, and they found only M. galloprovincialis. 
Similarly, our sample from the same site also yielded 
only M. galloprovincialis (Table I), but more extensive 
sampling inside San Diego Bay in a small-craft marina 
revealed both species and their hybrids. 

Although earlier reports indicated that hybridization 
was evident in San Francisco Bay (McDonald and Koehn, 
1988; Sarver and Foltz, 1993) and relatively infrequent 
elsewhere (Rawson and Hilbish, 1995), this survey re¬ 
vealed hybrids in San Diego Bay, Elkhorn Slough, two 
areas within San Francisco Bay, and two sites in Washing¬ 
ton. San Francisco Bay and Elkhorn Slough are located 
near the southern boundary for continuous populations of 
M. trossulus and the northern boundary for continuous 
populations of M. galloprovincialis, so hybridization 
there is not surprising. Studies of mussels in San Fran¬ 
cisco Harbor (Rawson et al., 1996) and the present study 
(Table I) revealed that the proportions of pure species, F, 
hybrids, and backcrosses vary across the bay, with A/. 
galloprovincialis becoming more common deeper in the 
bay. The studies of mussels from Palo Alto clearly indi¬ 
cate that hybridization and backcrossing are common, and 
that the mixing appears to be persistent over time. Pres¬ 
ence of hybrids away from this contact zone around San 
Francisco requires further explanation. The San Diego 
Bay site is in a marina where small craft may have intro¬ 
duced adult or larval M. trossulus, initiating hybridization. 
M. galloprovincialis and hybrids at Whidbey Island may 
be progeny of cultured mussels in a nearby mariculture 
facility. A hybrid at Tatoosh Island is possibly an indica¬ 
tion of transport of hybrid larvae from Puget Sound. 


MiTILUS DISTRIBUTIONS 


193 


Sporadic episodes of transplantation by means of ships 
will cause the presence of Alytilus species and their rela¬ 
tive abundances to fluctuate over time. Although the 
adults are effectively sedentary, the pelagic larvae are in 
the water column for 3 weeks or more; depending on the 
currents, they may move distances in excess of 200 km. 
But in addition to their natural dispersal, mussels are 
moved about extensively on the hulls of and in the ballast 
water of ships. Freighters moving between continents can 
inoculate a bay with many millions of larvae, and thus 
the distributions and the relative numbers of species of 
Alytilus might change from year to year (Carlton and 
Geller, 1993). 

The absence of M. edulis in our samples, as in previous 
allozyme and mtDNA reports, remains enigmatic, for 
there are several reports that Al. edulis has been trans¬ 
planted to the west coast of North America. For example, 
hundreds of oyster epibionts, which include Al. edulis, 
have been accidentally released in Pacific sites (Carlton. 
1989). Furthermore, we suspect that individuals of Al. 
edulis are introduced regularly to western ports by cargo 
ships discharging ballast water (Carlton, 1989; Carlton 
and Geller, 1993). Al. edulis has subsequently been identi¬ 
fied (D. Heath, pers. comm.) from a sample of alien mus¬ 
sels (Heath et al., 1995) collected in the Georges Strait, 
north and west of Vancouver, near an area where Al. 
edulis is farmed. 

Transition zones 

In California, the transition between Al. galloprovin- 
cialis and Al. trossulus is near 40° or 41°N latitude, ap¬ 
proximately the Cape Mendocino region. Previous studies 
had suggested that the northern limit of pure Al. gallo - 
provincialis was at Point Conception or Morrow 1 Bay 
(Sarver and Loudenslager, 1991; Rawson and Hilbish, 
1995), but our study revealed pure Al. galloprovineialis 
in Tomales Bay, nearly 700 km further north. This transi¬ 
tion is far north of the biogeographic boundary at Point 
Conception, a recognized demarcation between cold-wa¬ 
ter and warm-water biotas (Valentine, 1966), where many 
molluscan species have a distributional limit. 

Which variables control the latitudinal distributions of 
Al. galloprovineialis and Al. trossulus? Temperature is 
generally believed to control distributional boundaries of 
marine organisms (Valentine, 1966), and Alytilus is no 
exception (Seed, 1976). Al. edulis is able to withstand 
freezing in ground ice at — 20°C for 6 to 8 months of the 
year (Williams, 1970), but the northern distribution for 
A 1. edulis may be limited by the temperatures that permit 
breeding (e.g., Stubbings, 1954; Barnes, 1957). Similarly, 
the southern limit of the range of Al. edulis in the western 
North Atlantic is determined by summer temperatures, 
which are lethal for Al. edulis on the Outer Banks of 
North Carolina (Wells and Gray, 1960). 


In both the eastern and western Pacific, the transition 
zones between Al. galloprovineialis and Al. trossulus oc¬ 
cur where the annual mean temperature is between 13° 
and I4°C. On both sides of the Pacific, the southernmost 
populations of pure Al. trossulus (sites 9 and 22) experi¬ 
ence mean July temperatures of 13.39° and 13.74°C, re¬ 
spectively. We hypothesize that summer temperatures 
warmer than 13°-14°C approach the upper thermal limit 
for Al. trossulus and reduce the viabilities of larvae and 
adults. Although historically Al. trossulus was established 
south of this latitude in California, it was probably re¬ 
stricted to cooler sites. Geller (unpubl. data) used mito¬ 
chondrial 16S rRNA gene sequences to identify Al. trossu¬ 
lus collected from Catalina Island in 1900 and from Mon¬ 
terey Bay in the 1870s. Both Catalina Island and 
Monterey Bay are cooler than the collection sites em¬ 
ployed in the present study in San Diego Bay. We do not 
propose that the temperature of 14°C is an absolute bar¬ 
rier, for we sampled Al. trossulus from San Diego Bay 
(mean annual temperature = 16.02°C, mean July tempera¬ 
ture = 17.69°C). However, the presence of Al. trossulus 
within San Diego Bay but not on the open coast, at the 
Scripps Pier, suggests that persistent reintroduction by 
means of intracoastal boating is necessary to maintain 
Al. trossulus in the relatively warm localities in southern 
California. If our hypothesis is correct, Al. trossulus 
would not persist in San Diego Bay without recurrent 
introductions by ships. 
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